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Abstract 

A three-component dynamic system with influence of pumping and nonlinear dis- 
sipation describing a quantum cavity electrodynamic device is studied. Different 
dynamical regimes are investigated in terms of divergent trajectories approaches 
and fractal statistics. It has been shown, that in such a system stable and unstable 
dissipative structures type of limit cycles can be formed with variation of pumping 
and nonlinear dissipation rate. Transitions to chaotic regime and the corresponding 
chaotic attractor are studied in details. 
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1 Introduction 



The most intriguing phenomenon in nonhnear dynamical systems theory is 
a transition from regular dynamics to irregular one. Despite a transition to 
irregular dynamics can be driven by stochastic sources introduction into corre- 
sponding evolution equations, a special interest from the theoretical viewpoint 
is an observation of such regime in deterministic systems. A famous work by 
E. Lorenz [1] shows a possibility of regular dynamical systems to exhibit chaotic 
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regime. After, a lot of works concerning this problem allowed to develop math- 
ematical tools to study chaotic regime in different kind of physical, chemical, 
biological and other systems, see Refs. [2,4,3,5], for example. 

Mathematically studying of Lorenz system was performed in [6,7], where au- 
thors introduced the geometric Lorenz flows. Invariance principle for geometric 
Lorenz attractors were obtained in [8]. Description of dynamics of the Lorenz 
system having an invariant algebraic surface presented in [9]. 

In physics, thermodynamics in particular, chaos is applied in the study of tur- 
bulence leading to the understanding of self-organizing systems and system 
states (equilibrium, near equilibrium, the edge of chaos, and chaos). The oc- 
currence of a chaotic regime is related to the interplay between the instability 
and nonlinearity. The instability is responsible for the exponential divergence 
of two nearby trajectories, the nonlinearity bounds trajectories within the fi- 
nite volume of the phase space of the system. The combination these two 
mechanisms gives rise to a high sensitivity of the system to the initial condi- 
tions. 

It is well known, that optoelectronical devices are the simplest systems which 
manifest chaotic regime under deterministic conditions. In practice, a control- 
ling of dynamical regimes in optoelectronical devices is very actual task, be- 
cause different types of radiation can be used for different areas of human activ- 
ity [11,10]. The problem to control chaotic regimes in laser dynamics attracts 
much more attention in last decades. Nowadays new mechanisms which lead 
to chaotic regime generation were found [12,13,14]. Moreover, a rapidly de- 
veloping approach to study the quantum cavity electrodynamic models shows 
a possibility of chaotic radiation and the chaos control [15,16]. Prom a laser 
physics it follows that self-organization processes with dissipativc structures 
formation or generation of chaotic signal can be stimulated by introducing 
additional nonlinear medium into the Fabry-Perot cavity [21,10,17,20,19,18]. 
Experimental studies of the chaotic regime in such systems allows to generalize 
well known theoretical models, identifying additional nonlinearity sources, in 
order to provide a more detailed theoretical analysis and to predict different 
types of laser dynamics [23,24,25,22]. 

Prom the theoretical viewpoint, a problem to find ability of the system to 
manifest chaotic regimes and the control of chaos is of a special interest in 
a nonlinear dynamical system theory. Moreover, due to chaotic dynamics is 
accompanied by a different type of attractors appearing, another interesting 
problem is to investigate properties of such attractors. On physical grounds, 
the related actual problem is to investigate a chaotic regime in a generalized 
model of a dynamical system, whereas an initial model does not exhibits it 
itself. 
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In this paper we solve above problems for a three-parametric dynamical sys- 
tem realized in optoclcctronical devices. We are aimed to study a possibility 
of the well known two-level laser system type of Lorenz-Haken model [11,17], 
which can not show chaotic regime in actual ranges of system parameters, but 
manifests it if an additional nonlinear medium is introduced into the cavity. 
In our investigation we use suppositions of absorptive optical bistability and 
consider three-parametric system in order to show that a transition to chaotic 
regime can be controlled by the pumping and a bistability parameter propor- 
tional to the atomic density. The last one is related to nonhnear dependence 
of a relaxation time of one of dynamical variables versus its value. It will be 
shown that varying in a system parameters one can arrive at stable and unsta- 
ble dissipative structures, period doubling bifurcation and a chaotic regime. 
A corresponding chaotic attractor is studied by different approaches. 

The paper is organized in the following manner. In Section 2 we present a 
model of our system where we introduce theoretical constructions to model 
an influence of a nonlinear dissipation. Section 3 is devoted to development 
of the analytical approach to study possible regimes. In Section 4 we apply 
the derived formalism to investigate a reconstruction of the dynamics of the 
system. Section 5 is devoted to considering the main properties of the chaotic 
behaviour of a system. Main results and perspectives are collected in the 
Conclusions. 



2 Model 



Let us consider a three parametric dynamical model, which is widely used to 
describe self-organization processes: 

X ^ —X + y, 

(jy = -y + xz, (1) 
ez — {r — z) — xy, 

where dot means derivative respect to the time t. This model can be simply 
derived from the well known chaotic Lorenz system [1] using following rela- 
tions: t' = at, X = X/^/b, y = Y/y/b, z = r — Z , e = a/b. Here X, Y, Z are 
dynamical variables of the Lorenz system; a, b and r are related constants; we 
drop the prime in the time t for convenience. The dynamical system (1) can 
be obtained in the lasers theory where the Maxwell-type equations for electro- 
magnetic field and a density matrix evolution equation are exploited. It leads 
to the system of the Maxwell-Bloch type that is reduced to the Lorenz-Haken 
model (1) for two level laser system (see, for example, Refs. [17,10,26]). In such 
a case x, y, z are addressed to a strength of the electric field, polarization and 
an inversion population of energy levels. Models type of Eq.(l) are exten- 
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sively used to describe synergetic transitions in complex systems [27], where 
above variables x, y, z acquire a physical meaning of an order parameter, a 
conjugated field and a control parameter; the quantity r relates to the pump 
intensity and measures influence of the environment. 

A naive consideration of dynamical regimes in the model Eq.(l) at r > and 
(T ~ e ~ 1 shows that despite the system has nonlinearities in both second and 
third equations its behaviour is well defined. In this paper we consider a more 
general case reduced to the absorptive optical bistability models [20,17]. Prom 
the physical viewpoint absorptive optical bistability is related to possibility of 
the additional medium in the Fabry-Pcrot cavity (phthalocyanine fluid [18], 
gases SFq, BaCl^ and CO2 [19,17]) to transit a radiation with high intensities 
only with absorption of week signals. Despite models of such a type were con- 
sidered in the one-parametric case and only a stationary picture of dynamical 
system behaviour is analyzed (see Ref.[28] and citations therein), we will fo- 
cus onto describing possible dynamical regimes in the three-parametric model. 
Mathematically, our model reads 

i = -X + y + Jk, f> 

(jy=-y + xz, 
^ ez — {r — z) — xy. 

Here the additional force is introduced to take into account influence of 
an additional medium; k is the so-called bistability parameter proportional to 

the atomic density. Formally, it means a dispersion of a relaxation time for the 
variable x. Indeed, if we rewrite the first equation in the form x = —x/T{x)+y, 
then the dispersive relaxation time is r(a;) = 1 — k./ {1 + k, + x"^). 

Our main goal is to study possible scenarios of a dynamical regimes recon- 
struction in such three-parametric system under influence of the pumping r 
and the nonlinear dissipation controlled by k,. We will investigate conditions 
of a chaotic regime appearance. The corresponding irregular dynamics well be 
studied in details. It will be shown that a chaotic regime can be controlled by 
pumping and properties of an additional medium. 



3 Stability analysis 

Let us consider conditions where all variables are commensurable, setting 
0", e ~ 1. A behaviour of the system we present in the phase space (x, y, z). At 
first, let us study steady states. Setting x = 0, ?/ = and i = 0, one can find 
fixed points described by the values Xq'\ y^^^ and Zq'\ where /c = 1, 2, 3. The 



1 + ^2' 



(2) 
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first one has coordinates 



(3) 



and exists always; two symmetrical points, denoted with superscripts 2 and 3 
respectively, have coordinates 



(2,3) (2,3) (2,3) 



yd 



) ^0 



iVr — K — 1, 



(4) 



r — K 



r — K 



and realized under the condition 



r > K 



(5) 



It means, that Tc = k + 1 is a critical point for the bifurcation of doubling. 

According to the stationary states behaviour under influence of pumping and 
nonlinear dissipation, let us perform a local linear stability analysis [29,30]. 
Within the framework of the standard Lyapunov exponents approach time 
depended solutions of the system (2) are assumed in the form u oc e^*, A = 
A + iuo^ where u = (,t, y.z), A controls the stability of the stationary points, u 
determines a frequency, as usual. Magnitudes for real and imaginary parts of A 
arc calculated according to the Jacobi matrix elements Mj,- = idf^'^^ / duj ) , 

i,j = x,y,z; Z*^*^ represents the right hand side of the corresponding dynam- 
ical equations in Eq.(2); the subscript relates to the stationary value. As a 
result, the Jacobi matrix takes the form 



1 + K- 



Xn 



+ A 

2/0 



-1 

1 + A 

Xo 



-Xo 

1 + 



(6) 



Let us define the stability of the fixed point (3). A solution of the eigenvalue 
problem yields that the three eigenvalues are real and negative: 



Ai = Ai = -1< 0, 

A2 = As = --K - 1 - -V«;2 + 4r < 0, 

1 1 

A3 = A3 = --K - 1 + -VkF+It < 0. 



(7) 



It means that if both k, and r satisfy the condition (5), then the phase space 
is characterized by only node (3). 



5 



On the other hand, at r > Tc, all three fixed points are realized in the phase 
space. As it follows from our analysis the fixed point (3) changes its stability 
and becomes a saddle, due to A3 > 0. 



The stability of other two points with coordinates given by Eq.(4) can be set 
from solutions of the cubic equation 



+ 
+ 



3 + K 



2 + K-r 



(r — K,)^ 
2 + K-r 



A^ 



r ^ 2 + K, — r 
h 2k.- 



2 + K + K 



r — K (r — k)'^ 
2 + K-r\'^' _^ 



A 



(8) 



K 



Analytically, we can define only a border of the domain of system parameters 
where periodic solutions are realized. To consider a possibility of limit cycles 
formation in the phase space let us assume A = ±ia; and insert it into Eq(8). 
It yields the equation for such a border in the form 



2k 



2 + K 



-2k^ 



■ — kY 

'2 + K 



6 + 



r 1 



2r — K 



\ {r — kY J 

+4(k - 1) + 3r (— 1) = 0. 

\r — K J 



(9) 



As it will be seen below inside this domain in the plane (r, k) only stable 
periodic solutions are realized. 



4 Analysis of dynamical regimes 



To describe the dynamical regimes of the system under consideration we cal- 
culate the phase diagram shown in the plane (r, in Fig.l. Here the dashed 
line r = K + 1 divides the plane (r, k) into two domains denoted as 1 and 11 
and relates to the critical values of r and k. Inside the domain 1 only node is 
realized. The corresponding phase trajectories are shown in Fig. 2a. Above the 
dashed line a set of different dynamical regimes can be observed. Inside the 
domain II the node point transforms into a saddle point and two additional 
stable focuses appear (see Fig. 2b), the corresponding Lyapunov exponents has 
negative real parts, and equivalent imaginary ones. In the domain bounded by 
the solid line stable limit cycles are formed. This line is obtained as a solution 
of Eq.(9) with assumption that dynamical solutions of the system Eq.(2) are 
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periodic, i.e. A = ±ia;. When we move from the domain II to the domain III 
above focuses loss their stability and. as result, stable limit cycles are realized. 
All phase trajectories are attracted by these manifolds. The corresponding 
limit cycles with unstable focuses inside of them are shown in Fig. 2c. The 
domain IV bounded by solid, dotted and dash-dotted lines corresponds to the 
case where two stable and the corresponding two unstable limit cycles appear 
(see Fig. 2d). Such a transformation of dynamical system behaviour (transition 
from domain III into domain IV) is a result of double Hopf 's bifurcation. Here 
unstable focuses, shown in Fig. 2c, loss their unstabihty that leads to appearing 
unstable limit cycles inside stable ones (see Fig. 2d). When we cross the dotted 
line (transition from the domain II into domain IV) the stability of focuses 
shown in Fig. 2b is not changed, but two stable limit cycles and two unstable 
limit cycles appear, as Fig. 2d shows. A scenario of such a bifurcation is as 
follows. With an increase in one of two system parameters r or k in a point 
from the dotted line a new manifold type of semi-stable limit cycle is formed 
[2]. When we move to the domain IV this semi-stable limit cycle decomposes 
into the both outer stable and inner unstable limit cycles. Moving in opposite 
direction (IV — > II), one of the well known type of Hoph bifurcation occurs: 
an annihilation of both stable and unstable limit cycles. 

Let us consider a change of the system behaviour when we move from the 
domain III to the domain V bounded by solid and dash-dotted lines. If the 
dash-dotted line is crossed, then doubling period bifurcation occurs. Here two 
stable hmit cycles with unstable focuses, shown in Fig.2c, are transformed 
into the one limit cycle (Fig.2e), whose period equals two periods of the one 
cycle, shown in Fig. 2c. The behaviour of the system inside the domain V is 
the same as shown in Fig.2e. If we move from the domain V to the domain 
VI we get the period doubling scenario of chaos formation (see Fig.2f). If 
the hue that divides both domains V and VI is crossed, then the stability of 
focuses shown in Fig.2e is changed. It leads to a formation of two unstable 




O-' — 1 — i — 1 — I — ■ — I — > — i i ] — i — 

S 10 15 20 2S i!f 

Fig. 1. Phase diagram of the system (2) 
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Fig. 2. Phase portraits of the system (2) at : a) r = 
K = 5.0; c) r = 15.0, k = 10.0; d) r = 20.0, k = 
f) r = 35.0, K = 12.0; g) r = 39.0, k = 11.5 



10.0, K = 15.0; b) r = 15.0, 
9.0; e) r = 25.0, k = 15.0; 



limit cycles and the additional doubling period bifurcation. If we plunge into 
the dark zone the phase trajectories become irregular (see Fig.2g) due to the 
period doubling bifurcations. In such a case we obtain the chaotic attractor 
bounded by outside surface. Properties of such an attractor will be studied 
below. If we move from the domain IV to the domain VI the standard doubling 
period bifurcation occurs (see transition from Fig.2d to Fig.2f). The transition 
from the domain II into the domain VI is described by the following scenario. 
In the line that separates both domains II and VI a complex situation is 
realized. Here transverse intersections of both stable and unstable manifolds 
of the hyperbolic stationary point (saddle) occur. This indicates a homoclinic 
structure and stochastic layer appearance. In such a case the phase space is 
characterized by two unstable limit cycles and irregular behavior of trajectories 
outside of them (see Fig.2g). 

To elucidate that we deal with irregular behaviour let us use a quantity h 
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Fig. 3. The time evolution of a distance 6{t) between two nearest trajectories in the 
attractor Fig.2g at 5{to) = W'^VS 

known as Kolmogorov-Sinai entropy defined as follows 



lim 

5(ro) 



In 



r 



This quantity allows to observe behaviour of two trajectories 1 and 2 originat- 
ing at t = To from two starting points Ui(ro) and U2{tq). It is well known that 
if /i > then the system manifests chaotic motion. The corresponding calcula- 
tions of the separation 5(r) = \u2{t) — ui{t) \ between these trajectories versus 
time T when the trajectories are in attractor is shown in Fig.3. Defining the 
average inclination angle for the fluctuating curve, one can define Kolmogorov- 
Sinai entropy h = 0.1749. According to the definition of h the characteristic 
mixing time for the system tmix = h^^ = 5.7175 in the time units. At t ^ tmix 
the behaviour of a system is well predicted, whereas at t tmix only proba- 
bilistic description of the system dynamics is justified. From Fig.3 it follows 
that our attractor has a bounded size due to that distance 6{t) saturates at 
large time r. 

Despite the value of h might also depend on the choice of initial points of the 
trajectories it gives not only the qualitative, but also a quantitative specifi- 
cation of the dynamical regime. The more general criterion which takes into 
account the direction of the vector of the initial shift S{to) is the maximal 
(global) Lyapunov exponent defined in the form 

6u{t) 



A„ = A(5(ro)) = lim-ln 

I — >CiO ± 



where an upper limit and the norm \\u\\ are used. The quantity A^ allows 
to consider a divergence of any two initially close trajectories in phase space 
starting from points u{t) and u{t'). It follows, that the divergence of such 
trajectories is given by the dependence 5u{t) = Su{To)e^'^'^ . It is seen, that 
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r 



Fig. 4. Dependence of the maximal Lyapunov exponent vs. pumping r at k = 10.0 

at Am < trajectories of dynamical system converge and hence we get the 
regular dynamics, where the trajectories moves toward fixed points in the 
phase space. In the limit case A^, = one gets the limit cycle in the phase 
space. If Am > trajectories diverge and the system shows the chaotic regime. 
To determine the maximal Lyapunov exponent we use G.Benettin's algorithm 
[31]. The corresponding spectrum of maximal Lyapunov exponents A^ and 
the domain of its magnitudes are shown in Fig.l. 

To prove that above irregular behaviour of phase trajectories can be under- 
stood as chaos let us describe in details a behaviour of the maximal Lyapunov 
exponent A^ versus the pump intensity r at fixed k. The corresponding graph 
is shown in Fig.4. It is seen that in the domain of the pump intensity r G [0, r^] 
one gets A^ < 0. It means that in the phase space only node point is real- 
ized. If one chooses r G (r^,rj], then two stable focuses divided by the saddle 
point are realized with A^ < 0. In the domain r G {t'^iT'^] stable focuses 
transform into stable limit circles. One needs to note that at < r < r'^ 
the doubling period bifurcation takes place. At the such bifurcation point the 
theoretical value for the maximal Lyapunov exponent is A^ = 0. However, in 
computer simulations A^ 7^ 0, here a strong nonlinearity in the vicinity of the 
saddle point acts crucially and hence the obtained values A^ > should be 
understood as an artifact of computer simulations 0- The critical value rj? of 
the period doubling bifurcation in our system is obtained when the maximal 
Lyapunov exponent takes marginal value inside the above domain. In the in- 
terval r G {r^d -,^1.1^ we get A^ > that corresponds to divergence of phase 
trajectories and leads to chaotic regime formation. At r > rf^^ the phase space 
is characterized by stable focuses. 



To prove this we compare a similar doubling period bifurcation in Roessler system 
and have found that abruptly changes it values in the vicinity of such bifurcation 
point. 



10 




Fig. 5. Time dependence x{t) at r = 39.0, k = 11.5 from Fig.2g 




-0.2 J 



Fig. 6. Auto-correlation function of process x(t) at r = 39.0, k = 11.5. The corre- 
sponding dependence at small time is presented in the insertion 

5 Properties of the chaotic regime 

Let us consider statistical properties of irregular behaviour of the system, 
which corresponds to the phase portrait in Fig.2g in details. It is well known 
that a chaotic regime is characterized by fast decay of an auto-correlation 
function and continuous character of frequency spectrum. To show that our 
attractor acquires such properties let us calculate an auto-correlation function 
C(r) = (2/T) jQ^'^x{t)x{t + T)dt of the process x{t) shown in Fig.5. The auto- 
correlation function is presented in Fig. 6. It is seen, that auto-correlation 
function C(r) has fast falling down character as in chaotic systems. It means 
that the system losses its memory very quickly. Another criterion to show the 
chaotic regime is a continuous frequency spectrum. To obtain it let us use the 
Fourier transformation of the correlation function that gives power spectrum 
function S{uj) = yj^/ir /q^ C(r)e*'^'^dr. The corresponding frequency spectrum 
is shown in Fig. 7. It is seen, that spectral density function has a continuous 
character. Such a picture is typical for chaotic regimes [4]. 
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Fig. 8. The one-dimensional law Zn+i = f{zn) for the system (2) at r = 39.0, 
K = 11.5 

Finally, considering chaotic systems it is interesting to find a one-dimensional 
law that allows to represent genuine chaotic behaviour. To this end we cal- 
culate the first return map Zn+i = f{zn) by fixing the consecutive values of 
the maxima of the function z{t). Mathematically these maximal values are z- 
coordinates of points in the Poincare cross-section of the surface {r—z)—xy = 
(i = due to the third equation in the Eq.(2)). The corresponding one- 
dimensional law Zn+i = f{zn) is showu in Fig. 8. 

It is principally that points -Zn+i) he on the one-dimensional curve with 
an acute peak with a good accuracy. From the figure Fig. 8 it follows that if 
the initial value zi < z^, then the chaotic character of dependence Zn{n) is 
observed. The corresponding dependence z„ versus number of returns is shown 
in Fig. 9a, where chaotic character is visually seen. In the opposite case, when 
-2i > trajectories will converge to the fixed point z'^ as Fig. 9b shows. 

Other important physical quantity characterizing classical dynamic stochas- 
ticity is the fractal dimension of the attractor. For defining of fractal dimension 
we will use algorithm of Grassberger-Procaccia [32]. Let we have set of state 
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34.5- 



jj-j^ 3 6 9 n 

Fig. 9. Dependence Zn{n) for the different initial values zi with 
r = 39.0, K = 11.5: (a) Zi < z^, zi = 30.0; (b) zi > z^, zi = 36.0 



34.95 at 



vectors Ui, i = 1,2 ... ,N corresponding to the successive steps of numeric 
integration. In our case Ui is the complete set of variables (2) with the values 
corresponding to the moments of time t = ti. Then we can use numeric data 
for estimation of the following expression 



N 



C{e) = lim 



r,^^ N(N - 1) ,j 



where ^(v^) is a step function. According to the Grassberger-Procaccia method, 
fractal (correlation) dimension of the attractor may be defined as 



Dr — lim 



log C{s) 
log{e) 



An expected dependence of C{£) is e^". So, the corresponding plot in double 
logarithmic scale must be a line with angular coefficient Dc. Results of numeric 
calculations are presented in Fig. 10. A solid line corresponds to least-squares 
approximation of the results of data processing with ~ 1.84 ± 0.01. Fol- 
lowing Ruelle and Takens [33] one can say that our attractor can be identified 
as a strange chaotic attractor. 
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C(e) 




0.1 1 10 

e 



Fig. 10. Dependence of C(e) on the values of e, plotted by numerical integration of 
the set of equations (2) at r = 39.0, k = 11.5 

6 Conclusions 



We have studied different types of dynamical regimes for the three-component 
system with influence of pumping and nonlinear dissipation. It was shown, 
that system under consideration is related to absorptive optical bistability 
systems. Using Lyapunov exponents approach we have shown, that varying 
in system parameters different types of of stable dynamic behaviour of the 
system are realized. It was found, that controlling the pumping and/or prop- 
erties of additional nonlinear force, related to absorptive medium in optical 
systems, unstable and different types of stable dissipative structures can be 
formed. Considering dynamics in the three-dimensional phase space a set of 
bifurcations from the limit cycle is studied in details. 

It was shown, that in the system under consideration a strange chaotic at- 
tractor is realized only if additional nonlinear medium is introduced into the 
optical cavity. Properties of such a chaotic attractor were investigated using 
Kolmogorov- Sinai entropy, maximal Lyapunov exponents map, spectrum of 
the attractor, one- dimensional first return map, and fractal analysis. It was 
shown, that two neighbour trajectories in the attractor are divergent with 
positive Kolmogorov-Sinai entropy. The corresponding frequency spectrum is 
continuous, the temporal auto-corrclation function decays fast to zero. From 
the first return map it follows, that there is a unique initial point exist that 
separates regular and chaotic behaviour of the system. At last, a fractal (cor- 
relation) dimension of the attractor is Dc ~ 1.84±0.01. Prom above criteria on 
can conclude that obtained attractive manifold can be considered as a strange 
attractor as follows from Ruelle and Takens work [33]. Obtained chaotic regime 
can be controlled by nonlinear dissipation rate and pumping. Varying such pa- 
rameters we can obtain a stable focus and stable limit cycle from the chaotic 
regime. 
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Our results arc in good correspondence with well known results in systems with 
absorptive optical bistability and related to real experimental data for lasers of 
"C" -class [10,21]. It can be applied to study chaotic regimes in different kinds 
of dynamical system that can represent self-organization processes in models 
self-consistently described by hydrodynamic amplitude, conjugated field and 
control parameter in the form presented by the Lorenz-Haken model. 
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